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Abstract. In this contribution I show the physical limitations to the applica- 
tion of standard synthesis models from to 10 Gyr. I also present the multi- 
wavelength spectrum (from 7-rays to radio) of young star forming regions based 
on evolutionary population synthesis models that takes into account the stochas- 
ticity of the star formation proces. I show how the correlation of different ob- 
servables can allow to establish a better understanding of the star formation 
processes. Finally I compare the spectral energy distribution of AGNs and star 
forming regions and I show the possible connection of Star Forming Galaxies 
with Seyfert 2 AGNs. 



1. Introduction: Intrinsic limits of synthesis models 

Since the work of Tinsley & Gunn (1976), evolutionary synthesis models have been 
extensively used to obtain the integrated properties of systems where stars are formed. 
Models are based on the convolution of isochrones with the Initial Mass Function (IMF) 
and the Star Formation History. For the case of an instantaneous burst of star formation 
(simple stellar population, SSP), the mean luminosity in a given band and a given age, 
fJ-L(t), results from the weighted sum of the number of stars with initial mass to^, Wi, 
(given by the IMF), and the individual luminosities li{mi,t) (given by the isochrone). 
If the sum of the Wi values is normalized to 1 Mq transformed into stars from the 
onset of the burst (as usual), the resulting luminosity will also be normalized. The total 
luminosity of a cluster is then, directly proportional to the initial mass transformed into 
stars, M: L^'"^(t) = x ^i(t). 

Such a modeling has some intrinsic constraints: 

• The total luminosity of the cluster modeled, L'^'"*(t), must be larger than the 
individual contribution of any of the stars included in the model, and, in particular, 
larger than the most luminous star, /™"^ (771, , t) . This statement defines a natural 
theoretical limit, that is not always considered when the models are applied to real 
observations. The modeling of clusters with masses below this limit can only be 
performed by codes that include sampling effects, either analytically or via Monte 
Carlo simulations (see, e.g. Cervifio et al. 2002). 

• There is also an uncertainty coming from the very nature of the IMF. Most the- 
oretical models assume that the IMF is completely populated, but Nature does 
not follow such rule. In other words, any modeling that assumes a completely 
populated IMF will be correct only under the asymptotic assumption of an infi- 
nite number of stars. Otherwise, the modeling only yields a mean value of the 
observed quantities. In fact, Monte Carlo simulations for clusters at different ages 
show that the mean values of the synthesized quantities depend on the amount 
of stars used for the simulations (see Santos & Frogel 1997, Cervino, Luridiana 
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& Castander 2000 or Bruzual 2002 as examples). The main difference between 
Monte Carlo and standard simulations is that Monte Carlo simulations always use 
an integer number of stars. On the other hand, by construction, standard (ana- 
lytical) simulations always assume that most of the mass values are represented 
by a fractional number of stars. 

Note that such uncertainty is relevant to any modeling that makes use of the IMF. 
An example of its influence in galactic chemical evolution models can be found in 
Cervifio & MoUa (2002). 



safe results 
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Figure 1. M ranges where different situations may occur when the results 
of synthesis models are compared with observed data as a function of time and 
metalhcity: Z=0.019 (solid line), Z=0.008 (dashed hne) and Z==0.0004 (point- 
dashed line). The upper line defines the value where there are, at least, 
10 post-main sequence stars (t < 10^ yr) or where i^o;* > 10 x /^°^(mi,t) 
(t > 10«). 

Then, before using a synthesis model, it is necessary to know the A4 value for which 
a variation of ±1 star is relevant to the results. This value can be estimated imposing 
that i^oT^ larger than 10 x lj^f^{mi, t). The resulting values for Ai for a Salpeter IMF 
in the mass range 0.08 - 120 Mq are shown in Fig. 1 using the results from Girardi et 
al. (2000)0. It yields a value of M > 10^ Mq, that assure that there are, at least, 10 
Post-main sequence stars {t < 10^ yr) or that Lg^"'* > 10 x /^°^(mj,i) {t > 10^). 



Available at 



ittp : //pleiadi . pd . astro .it/ 
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2. Statistical synthesis models 

The only way to deal with small systems is to include statistical effects in the synthesis 
models. This can be performed computing the variance of the observable, cr£(i), to- 
gether with the mean value. It can be done easily assuming that the Wi values follow a 
Poissonian distribution (they must always be positive and integer numbers, see Cerviiio 
et al. 2002 for details). However, it is more useful to show the variance in terms of an 
effective number of stars, ML{t) = ^i\{t)/(j\{t) which was defined by Buzzoni (1989). 
Both, variance and TV, scale linearly with M: the smaller the system, the smaller TV 
(and the larger the relative error). This formulation has the additional advantage of en- 
abling us to obtain the correlation coefficient between two observables, p{Li, L2). This 
coefficient is independent of TW and it only depends on the IMF slope, so it can also 
be used for the analysis of large samples (see the contribution from R. Terlevich on the 
subject of mega datasets). Finally, the use of TV also enables us to know when multi- 
or bi-modal distributions may be observed and when the results of synthesis models 
are biased in comparison with the observations (see Cervifio & Valls-Gabaud 2002 for 
details). 

In the following I am going to apply our synthesis models results to different astro- 
physical scenarii. Model outputs^ can be found in the www server 

ittp : //www . laef f . esa. es/users/mcs/SED 

2.1. Multiwavelength emission 

In Fig. 2 I show the multiwavelength spectrum of a 5.5 Myr old burst of star formation 
including the 90% confidence interval that arise from sampling fluctuations in the stellar 
population for clusters with a given M and a Salpeter IMF slope with mass limits 2-120 
Mq . It is interesting to note that the lower dispersion corresponds to the UV continuum, 
which turns out to be the most reliable age indicator (in absence of extinction effects). 

Additionally, the correlation coefficients obtained theoretically can allow to estab- 
lish which are the better observables to determine physical properties. Let me give you 
an example: The presence of Wolf-Rayet (WR) stars imply the presence of ionizing 
flux, hence, some degree of correlation is expected in the age-diagnostic diagrams of 
EW(H/3) vs. EW(WR). The age-dependent correlation coefficient defines the geometry 
and the way how the observational data must be interpreted in such a plot, in the sense 
that variations in one of the observables also determine how the other observable varies. 
Such variations are only dependent on the age and the IMF slope since the correlation 
coefficient is sampling independent. 

Finally, there remains the problem of how to apply synthesis models to individual 
small systems. In this case, the situation is more complicated. Taking again the example 
of WR stars, whereas the presence of such stars enable us to quote an age range where 
the star formation began, their absence cannot be used to quote any plausible age 
range. However, if synthesis models have a physical reality, different observables must 
be related to each other. Hence, the real solution will lie in the joint distribution of the 
probability distributions of the different observables. Unfortunately, this kind of study 
has not yet been performed. 

2.2. Starburts-AGN connection 

In the preceding items I have shown that synthesis models must be applied with cautions 
to small systems. However, the definition of small depends on the statistics of the 
observable. In Cervifio, Mas-Hesse & Kunth (2002) we have computed the soft X-ray 
emission from starburst galaxies. This emission is produced from the reprocessing of the 
kinetic energy produced by the star forming regions and from the Supernova Remnants. 



^Not all the model output are presented in the www server, please, ask us for any special request. 
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In this case, the dispersion arises mainly from the occurrence of Supernova (SN) events, 
that is intrinsically a poor statistical observable: the SN rate is around 10~^ SN yr"-'^ 
Mq^ (following a Salpeter IMF in the mass range 2 - 120 Mq, Cerviho & Mas-Hesse 
1994), and the X-ray life time of an event is around lO'' yr, which implies Ai larger 
than 10"' in order to have a mean value of 1 SN active in X-rays. In that work the X-ray 
emission of individual SN events was averaged over the SN life time (see Aretxaga et al. 
in this proceedings for the statistics and evolution of SN fluxes) , hence the quoted values 
are an average estimation with an additional dispersion that has not been considered. 
This additional dispersion may be relevant in the interpretation of the optical emission 
line spectrum in star forming galaxies (see also Rodrfguez-Gaspar & Tenorio-Tagle 1998, 
Silich et al 2001, and Stasinska & Izotov 2002) 

Besides the limitations, the approach is enough to test the relevance of circumnu- 
clear star forming regions in the energy budget of Active Galactic nuclei ( AGN) . In Fig. 
3, left panel, I show the soft-X ray/UV ratio predicted by the synthesis models and 
their comparison with the observational data of Star forming galaxies. In all the cases, 
some amount of kinetic energy must be reprocessed in X-ray in order to explain the 
observations. The comparison with Seyfert galaxies is shown in the right panel of the 
figure. Star forming regions are able to explain the X-ray emission of some Seyfert 2 
galaxies, but their emission is not enough to explain the X-ray emission in Seyfert 1. In 
this approach, the putative black hole in some Seyfert 2 galaxies will only be relevant 
in the hard X-ray domain. 
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Figure 3. Left: Ratio L^(0.1 - 3.5 keV) over z/Li^(1450 A) computed for a 
20% reprocessing values of mechanical energy into soft X rays and different 
metallicities. In grey, the 90% CL for a cluster where 10^ M0 have been trans- 
formed into stars for two extreme metallicities and reprocessing values. Data 
points from star forming galaxies are marked as triangles. Right: Histograms 
for the ratio X-ray/UV for different types of emission line galaxies. The ver- 
tical solid line corresponds to the upper value predicted by our evolutionary 
synthesis models. 
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Discussion 



Dottori: How circumnuclear are the circumnuclear star forming regions you are talking 
about? 

Cervifio: The data I have presented are based on lUE observations, so they are in a 
region of 10 x 20 arc sec around the center. 

Cid Fernandes: You have rightfully warned us that evolutionary synthesis predictions 
come with a statistical uncertainty. Have you compared your predictions for the vari- 
ances with the observed dispersion in properties of star-forming regions? 

Cerviiio: Not personally. At this moment I am still working on the statistical impli- 
cations. However, other authors (Bresolin & Kennicutt 2002) have applied the models 
results available on the Web server to high-metallicity star-forming regions and the the- 
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oretical dispersion looks to be consistent with the observed data. 

Garcia Vargas: It is true that taking into account all the probabilistic effects on clus- 
ter modeling is important to understand the dispersion, or degree of correlation when 
plotting observed data of GEHRs. However I think that this approach should be sys- 
tematically applied to large GEHR where there is evidence of several small clusters 
(< 10^ from the region (as revealed from HST) 

Cervino: Yes, indeed the effect of sampling is more dramatic for low mass clusters 
{A4 < 10^ M0), but it will be also present in more massive systems. 



